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X-Ray Observations of the W51 Complex with Suzaku 
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Abstract 

We present a detailed analysis of the X-ray emission from the middle-aged supernova remnant W51C 
and star-forming region W51B with Suzaku. The soft X-ray emission from W51C is well represented by 
an optically thin thermal plasma in the non-equilibrium ionization state with a temperature of ~0.7 keV. 
The elemental abundance of Mg is significantly higher than the solar value. We find no significant feature 
of an over- ionized plasma in W51C. The hard X-ray emission is spatially coincident with the molecular 
clouds associated with W51B, overlapping with W51C. The spectrum is represented by an optically thin 
thermal plasma with a temperature of ~5 keV or a powerlaw model with a photon index of ~2.2. The 
emission probably has diffuse nature since its luminosity of 1x10'^"' erg s^^ in the 0.5-10 keV band cannot 
be explained by the emission from point sources in this region. We discuss the possibility that the hard 
X-ray emission comes from stellar winds of OB stars in W51B or accelerated particles in W51C. 

Key words: ISM: individual (W51C) — ISM: supernova remnants — ISM: cosmic rays — ISM: Hll 
regions — X-rays: ISM 



1. Introduction 

Galactic cosmic rays (mostly protons) are widely be- 
lieved to be accelerated at the shock of supernova rem- 
nants (SNRs) through the diffusive shock acceleration 
mechanism (Reynolds 2008, and references therein). TeV 
7-ray observations of SNRs gave us firm evidences of accel- 
eration of TeV particles (e.g., Aharonian et al. 2004). In 
addition, the signatures of proton acceleration were ob- 
tained by GeV 7-ray observations of middle-aged SNRs 
interacting with molecular clouds (e.g., Abdo et al. 2009). 
On the other hand, how the acceleration process evolves 
in SNRs is currently not well understood although it is 
crucial in understanding the acceleration mechanism in 
SNRs. 

Theoretically, one possible scenario of the evolution of 
acceleration in SNRs is that eflticient acceleration, acceler- 
ating a large amount of particles up to high energy, occurs 
at the free expansion phase and it ends at the early Sedov 
phase. The cosmic-ray streaming instability becomes less 
efficient as the shock velocity decreases with time, and 
nonlinear wave interactions reduce the level of turbulence 
at the late Sedov stage (Voelk et al. 1988). An alternative 
scenario is that moderate acceleration lasts until the end 
of the Sedov phase. In this case, tens-of-TeV particles can 
exist in middle-aged SNRs (Sturner et al. 1997; Nakamura 
et al. 2012). Thus it is important to investigate nonther- 
mal particles in middle-aged SNRs. 

The spectrum of high energy electrons traced by syn- 
chrotron X-rays could have an important clue of the evo- 



lution of the acceleration mechanism since its cooling time 
is fairly shorter than that of protons and the electron 
spectrum can be strongly affected by the cooling process. 
Possible synchrotron X-ray emission from a middle-aged 
SNR was detected in W44 (Uchida et al. 2013). However, 
the number of the samples is limited because emission 
from the thermal plasma of middle-aged SNRs dominates 
in the X-ray band and nonthermal emission is buried in it. 
Thus, accurate measurement of the thermal component is 
required to detect the nonthermal emission, in addition to 
observations with a highly sensitive detector in the hard 
X-ray band. 

Moreover, investigations of SNR plasmas can also 
constrain the time evolution of particle acceleration. 
Suprathermal particles injected into the cosmic-ray accel- 
eration can effectively ionize ions in plasma and generate 
an over-ionized plasma which has a higher ionization tem- 
perature (Tz) than the electron one (To: Kato & Masai 
1992). In big solar flares, the ionization rate of plasma 
increases as the flux of hard X-rays from bremsstrahlung 
emission increases (Tanaka 1986). The degree of over- 
ionization depends on the elapsed time of the ionization 
and the density of suprathermal electrons. Thus, we can 
constrain the density of accelerated cosmic rays and the 
time when the acceleration occurred in the SNR via the 
diagnosis of thermal plasma with X-ray. 

W51C (G49.2-0.7) is an SNR with the age of 
^3x10'^ yr derived from a Sedov analysis with the ob- 
served X-ray temperature and the angular radius (Koo et 
al. 1995). The extended 7-ray emission has been detected 



2 



Yoshitaka Hanabata et al. 



by the Fermi satellite, the H.E.S.S. and the MAGIC tele- 
scopes (Abdo et al. 2009; Fiasson et al. 2009; Aleksic et 
al. 2012). The 7 rays arc naturally explained by Tr^-decay 
through the interaction of protons accelerated in the SNR 
with dense gas (Abdo et al. 2009; Aleksic et al. 2012). 
Thus, W51C is a good SNR to study cosmic-ray accel- 
eration during the Sedov phase. This SNR is located at 
the tangential point of the Sagittarius Arm where the dis- 
tance is estimated to be 6 kpc (Koo et al. 1995). The 
remnant appears as an incomplete shell of ~30' diame- 
ter in the radio continuum (Subrahmanyan & Goss 1995). 
The western part of the remnant is superposed upon the 
massive star-forming region W51B (Koo 1999; Kumar et 
al. 2004). Koo & Moon (1997) found that W51G inter- 
acts with high-velocity molecular stream and Carpenter 
& Sanders (1998) revealed that it most likely belongs to 
the cloud associated with Hll regions of W51B by the 
CO observation. The detections of two OH (1720 MHz) 
maser spots give us another evidence of interaction be- 
tween W51C and a part of molecular cloud associated with 
W51B (Green et al. 1997). 

X-ray observations of W51C were carried out by 
Einstein, ROSAT, ASCA and Chandra (Seward 1990; Koo 
et al. 1995; Koo et al. 2002; Koo et al. 2005). The soft X- 
ray emission from the SNR was represented by the plasma 
emission in collisional ionization equilibrium (CIE) with 
the temperature of 0.3-0.5 keV. No hint of the over-ionized 
plasma has been reported from W51C. In the hard X-ray 
band, ASCA and Chandra detected two pulsar wind neb- 
ula (PWN) candidates and some H II regions in addition to 
the extended emission overlapped with molecular clouds 
of W51B (Koo et al. 2002; Koo et al. 2005). However, 
because of the low sensitivity of ASCA and Chandra for 
hard X-rays, the spectral analysis for the extended emis- 
sion was not performed and its origin is unclear. 

The aim of this paper is to search for the over-ionized 
plasma and reveal the origin of the extended hard X- 
rays for studying the evolution of acceleration process in 
middle-aged SNRs. This study was done with the X-ray 
Imaging Spectrometer (XIS: Koyamaet al. 2007) onboard 
Suzaku (Mitsuda et al. 2007) with the superior energy res- 
olution and the large effective area even in the hard X-ray 
band. 

2. Observation and Data Reduction 

We carried out two observations around the W51C re- 
gion. Locations of two pointing fields and an observation 
log are shown in figure 1 and table 1, respectively. 

The XIS consists of four X-ray CCD cameras, 
each mounted on the focal plane of the X-ray 
Telescope (XRT: Scrlemitsos et al. 2007). Three sensors 
employ Front-IUuminated (FI) CCDs (XIS 0, 2 and 3) 
while the other employs a Back-Illuminated (BI) CCD 
(XIS 1). The entire region of XIS 2 and an edge of XIS 
have not been functional due to the anomalies in 2006 
November and in 2009 June, respectively. The XRT has 
a point spread function of ^2'. The total effective areas 
is about 400 cm^ at 8 keV. 




Fig. 1. The ROSAT image (gray-scale; 0.7-2.0 keV) of 
W51C region (Koo et al. 1995). The thick and thin solid 
boxes indicate the position of XIS and Chandra field of views, 
respectively. The blue contours show the VLA 1.4 GHz radio 
image (Koo & Moon 1997) with a level of 0.015, 0.06, 0.2, 0.6 
and 1.5 Jy boam~^. 

The XIS was operated in the normal clocking full- 
window mode, and the data of the two editing modes, 
3x3 and 5x5, were combined. To restore the radiation- 
induced degradation in the energy gain and resolution, 
the spaced-row charge injection technique (Bautz et al. 
2004) was applied with the makcpi files version 20110708 
provided by the XIS team (Uchiyama et al. 2009). In the 
analysis, we employed the data with a process version 2.5 
for W51C_W and W51C_E and a version 2.4 for the back- 
ground data, respectively. The data were analyzed with 
the HEAsoft version 6.11. We removed hot and fiicker- 
ing pixels, and discarded events with a condition during 
the South Atlantic Anomaly passages. Earth night-time 
elevation angle below 5°, and day-time angle below 20°. 

We also used Chandra archival data of observations 
around W51C to evaluate the contribution of point 
sources to the extended hard X-rays. The observation 
log is also shown in table 1. 

3. Analysis 

3.1. Imaging Analysis 

Figure 2 shows the XIS mosaic images around W51C 
in the energy bands of 0.5-2.5 keV and 3.5-8.0 kcV using 
XIS 1 and 3. The non- X-ray background (NXB) created 
by xisnxbgen (Tawa et al. 2008) was subtracted and then 
vignetting and exposure corrections were performed. 

The W51C_E observation covers the central bright 
region of the X-ray emission from W51C, while the 
W51C_W covers the large portion of the star-forming re- 
gion W51B (see figure 1). The 0.5-2.5 keV image shows 
an extended emission along the east-center direction in 
the field, which has two bright portions located at south- 
east and central north (figure 2a; Reg 1 and Reg 2). 
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Table 1. Log of Suzaku and Chandra observations of W51C region. 



Instrument 


Field name 


Sequence no. 


Aim point* 


Start date 


Effective 








m 


(UTC) 


exposure 


Suzaku 


W51C.W 


504066010 


(49.°11,-0.°31) 


2010/03/28 


44.1 ks 




W51C_E 


504067010 


(49.°ll,-0.°54) 


2010/03/30 


43.7 ks 




BCD 


504044010 


(65.°87,-0.°32) 


2009/05/23 


38.7 ks 


Chandra 


W51B 


500318 


(48.°96,-0.°31) 


2003/06/03 


29.9 ks 




W51_N0RTH 


500319 


(49.°17,-0.°35) 


2002/12/08 


11.8 ks 



* In Galactic coordinates. 




0.5 1.0 1.5 2.0 2.5 0.2 0.4 0.6 0.8 1.0 1.2 

(xlO"^ counts s"^ pixel"^) (xlO'^ counts s"^ pixel'^) 



Fig. 2. Suzaku XIS counts map of W51C region in the (a) 0.5-2.5 keV and (b) 3.5-8.0 keV bands. Gray contours show ^^CO 
(J = 2—1) line intensity (Bieging et al. 2010) at 10, 20 and 40% levels, for the velocity range from 60 to 70 km s^^. Red contours 
indicate the intensity distribution of the high velocity Hi gas at 40 and 80% levels, integrated in the velocity range from 85 to 
144 km s~^. Black ellipses indicate the regions used for the spectral analyses. Dashed regions indicate the excluded regions for the 
analysis, where two blue ones, a north large green one, and the other green ones are PWN candidates, calibration source, and Hll 
region of G49. 0—0.3 and G49.2— 0.3, respectively. Two green crosses indicate the Hll region of G49. 1—0.4 and G49. 10—0.27. 



These structures are consistent with the images obtained 
by ROSAT (Koo et al. 1995) and ASCA 0.7-2.5 keV 
band (Koo et al. 2002). 

W51C interacts with some part of the molecular cloud 
to which W51B most likely belongs (Koo & Moon 1997; 
Green et al. 1997). To investigate the correlation be- 
tween the X-rays and the molecular clouds, we overlaid the 
^^CO(J = 2-1) line intensity map (Bieging et al. 2010) for 
the velocity range from 60 to 70 km s~^ on the XIS images. 
In the soft X-ray image, the surface brightness drops down 
toward the western edge of the clouds. Koo et al. (1995) 
and Koo et al. (2002) concluded that it was caused by 
the absorption of the intervening clouds. The 3.5-8.0 keV 
image shows an extended (~8'xl2') emission, which over- 
laps with the molecular clouds at the western part of the 
SNR (figure 2b; Reg 3). There are also two PWN can- 
didates and HII regions G49.0-0.3 and G49.2-0.3 in the 
hard X-ray image. In this paper, we focus on the emission 
from W51C and the extended hard X-rays. 

We also compared the X-ray emission with the 



shocked Hi clouds with a velocity range from 85 to 
144 km s-i (Koo & Moon 1997). The clouds are lo- 
cated near the western boundary of the soft X-rays, while 
the correlation between the hard X-ray emission and the 
clouds is not clear. 

We first investigate the spectra of bright regions in the 
soft X-ray band (subsection 3.2); Reg 1 and Reg 2, which 
could be the parts of the SNR. Then we inspect the spec- 
trum of the hard X-ray emission in Reg 3 (subsection 3.3). 

In the spectral analyses, we used the SPEX soft- 
ware (Kaastra et al. 1996) version 2.02.04. The redistri- 
bution matrix and the auxiliary response functions were 
generated by xisrmfgen and xissimarf gen (Ishisaki et 
al. 2007), respectively, with the rmfparam files released on 
20110804. The errors are quoted at the 90% confidence 
level unless otherwise mentioned. 
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Fig. 3. Background-subtracted spectra with the best-fit models. Each panel shows the result of Reg 1 for the CIE (a) and the 
NEI (b) model and Reg 2 for the CIE (c) and the NIE (d) model. Black and red data show the FI and BI spectra, respectively. 
For visibility, the FI spectrum is multiplied by 5. The vertical error bar of each data point represents the la error. The Gaussian 
component for FI and BI are shown with green solid and dashed lines, respectively. 



3.2. Spectral Analysis of Soft X-Ray Emission 

3.2.1. Background Estimation 

There is no observation of a blank sky nearby W51C. To 
estimate the background for the soft X-ray emission care- 
fully, we sought a Suzaku archival dataset satisfying the 
following criteria: (i) 40° < Z < 70°, -0.°50 <h< -0.°20, 
(ii) exposure > 40 ks, (iii) free from bright sources, and 
(iv) observed with normal clocking mode. As a result, we 
found that the background field (hereafter, BGD) listed 
in table 1 only meets the criteria. 

The background emission of XIS consists of three com- 
ponents, i.e. NXB, cosmic X-ray background (CXB) and 
Galactic ridge X-ray emission (GRXE: e.g., Koyama et al. 
1989). NXB is estimated by using xisnxbgen as described 
in subsection 3.1. Since CXB has the uniform spectral 
shape and intensity in all sky (Kushino et al. 2002), we can 
model the CXB spectrum in an arbitrary region. GRXE 
has strong positional dependencies. The longitudinal and 
latitudinal intensity profiles of the Galactic center X-ray 
emission and GRXE Fe XXV Ka (6.7 kcV) line were mod- 
eled by Uchiyama et al. (2013), which were obtained by 
the analysis of XIS data with |/| < 125° and \b\ < 5°. Using 



this model, wc can derive the GRXE intensity ratio be- 
tween two arbitrary regions. Note that we did not use 
the result of Kaneda et al. (1997) which was obtained 
by the study of the latitudinal dependence of GRXE in 
the Scutum arm region (28° < I < 29°). This is because 
they could not resolve the Pel Ka (6.4 keV), XXV Ka 
and XXVI Lya (6.9 keV) lines and represented the spec- 
trum of GRXE with one foreground emission (FE) plus 
one thermal plasma component, which is insufficient (see 
subsection 3.3.1). 

For the analysis of Reg 1 and Reg 2, we modeled 
the GRXE spectrum toward W51C by scahng that to- 
ward the BGD observation region. We first excluded 
the point sources and the region of 6 > —(f.20 from 
BGD because the analysis regions in W51C are lo- 
cated at b < — 0?20 (see figure 2). The solid angle of 
the remaining region was 221 arcmin^. Then we sub- 
tracted NXB and CXB to create the GRXE spectrum 
toward the BGD region. We assumed the CXB spec- 
trum to have a powerlaw with a photon index of 1.41 
and an intensity of 5.4x10"^^ erg s~^ arcmin"^ in 2- 
10 keV band (Kushino et al. 2002). The absorption 
of interstellar medium (ISM: Morrison & McCammon 



No. 



X-Ray Observations of W51 with Suzaku 

Table 2. The best-fit results* of Reg 1 and Reg 2. 



Parameter 



Reg 1 



Reg 2 



CIE 

iVn (10^2 cm 
VEM (10'" 
To (keV) 
Zug (solar) 
Zsi (solar) 
Zs (solar) 
xVd.o.f 



1.56 (1.61) 
7.68 (7.74) 
0.49 (0.48) 
1.91 (2.02) 
1.08 (1.14) 
1.46 (1.49) 
448/348 



+0.06 
-0.05 
+ 1.31 
-0.96 
+0.02 
-0.02 
+0.18 
-0.17 
+0.13 
-0.12 
+0.32 
-0.29 



2.01 (2.07) 
7.17 (7.47) 
0.56 (0.55) 
1.53 (1.61) 
0.50 (0.51) 
0.40 (0.40) 
200/160 



+0.09 
-0.08 
+ 1.41 
-1.01 
+0.02 
-0.03 
+0.20 
-0.20 
+0.09 
-0.09 
+0.19 
-0.18 



NEI 

Nh (10^2 
VEM (10'" 
Te (keV) 



~3\t 
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1.49 (1.51) t"o"ol 

4.14 (3.94) +0;^° 

0.70 (0.69) t'ot 

1.77 (1.89) t^o.ll 

1.73 (1.91) tVil 

0.89 (0.97) t°o.ll 

0.83 (0.88) toil 
383/347 



Tlet (10 
Zug (solar) 
Zsi (solar) 
Zs (solar) 

xVd.o.f 

* 'I'he values in the parentheses are obtained by increasing the 
GRXE spectrum by a factor of 1.25, for systematics. 
t Volume emission measure, VEM = J ncUYidV, where rio and 
riH are the electron and Hydrogen densities, respectively, and V 
is the emitting volume, assuming the distance of 6 kpc. 



1.97 (2.02) torn 
5.06 (5.12) % 
0.69 (0.68) t^ol 
2.97 (3.14) 
1.40 (1.50) toll 
0.48 (0.50) 
0.30 (0.30) toil 
168/159 



1983) for CXB is attributed to Galactic Hi medium of 
-0.8x10^2 cm-2 (Kalberla et al. 2005)^ and molecu- 
lar hydrogen of -0.2x10^2 cm-^ (Dame et al. 2001). 
Then, the overall absorption column density for CXB is 
^tptai ^ 7VHi + 2iVH2 = 1.2x10^2 cm-2. We assumed the 
solar abundance for the absorbing materials (Anders & 
Grevesse 1989). 

To derive a GRXE scaling factor of the W51G region 
against the BGD region, we used the GRXE intensity pro- 
file and then the scaling factors were obtained to be 1.76 
and 1.86 for Reg 1 and Reg 2, respectively. Since the 
profile of the GRXE intensity was estimated with an un- 
certainty of 25%, we take into account it as the systematic 
error for the spectral analysis. 
3.2.2. Reg 1 and Reg 2 

To investigate the nature of the soft X-ray emission, 
we extracted the spectra from Reg 1 and Reg 2. For 
background subtraction, we first subtracted NXB and 
CXB considering the ISM absorption with the A?^^°'^' of 
1.9x10^2 cm-2 and 2.8x10^2 cm-2 for Reg 1 and Reg 2, 
respectively. These column densities are determined in 
the same manner as subsection 3.2.1. Then the scaled 
GRXE spectrum prepared in subsection 3.2.1 was sub- 
tracted from each region after correcting the effective area. 
The background-subtracted spectra are shown in figure 3. 
We can see Ka emission lines from He-like (hereafter, 
Hea) neon (Ne). magnesium (Mg), silicon (Si) and sul- 
fur (S). 

For Reg 1, we first tried to fit the spectrum with an 

^ We used the Hi calculator available at 
(http: / /heasarc. nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi). 



optically thin thermal plasma in CIE with an ISM ab- 
sorption. The abundances of the noticeable elements Mg, 
Si and S {Zyig, Z^i and Zg) were free parameters. Those 
of Ne and Fe were fixed at the solar value because the Ne- 
K and Fe-L lines are overlapping and the heavy absorp- 
tion make it hard to determine them independently. The 
other elements were also fixed at the solar values. This 
model failed to reproduce the spectrum with a x^/d.o.f of 
593/349, leaving large residuals around 0.9 keV, 1.2 keV, 
and > 3.0 keV. The first residual around 0.9 keV corre- 
sponds to the Ne Hea. The second is probably due to the 
lack of a number of Fe L-shell lines emitted at the tran- 
sitions from highly excited states of n > 5 (Brickhouse et 
al. 2000) in the plasma model. Thus, a Gaussian model 
was artificially added at 1.2 keV. The fit was largely im- 
proved to a ^ jd.o.l of 448/348 although the residuals 
around 0.9 keV and > 3 keV still remain. The best-fit 
spectra and parameters for the one CIE model are given 
in figure 3a and table 2, respectively. The value of 
was obtained to be 1.9lio |^y solar, which is significantly 
higher than the solar value. 

We next examined a non-equilibrium ionization (NEI) 
model. In the NEI model, the ionization parameter n^t 
was an additional free parameter, where rio and t arc the 
electron density and the elapsed time after the shock heat- 
ing. The fit was significantly improved to a 'x^ jA.o.i of 
383/347 without a notable residual around 0.9 keV (fig- 
ure 3b and table 2). The values of and riat are obtained 
to be 0.70t[J;E!| kcV and 1.771^:42 xlO" cm-^ s, respec- 
tively. The obtained of 1.731q ^[4 is still significantly 
higher than the solar value. 
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We also checked whether the plasma is over-ionized 
since the excess around 0.9 keV and > 3 keV could be 
caused by the radiative recombination continua (RRCs). 
RRCs appear when the recombination process dominates 
over the ionization one in the over-ionized plasma (T^ > 
Tg). If RRCs generate the notable residuals, the excesses 
around 0.9 keV and > 3 keV would correspond to the K- 
shell recombinations of H-like O and Si ions, respectively. 
Hence, we set Zq and T^. to be free in addition to the above 
CIE model. As a result, we obtained Tz/Tc of 0.83i5]:n, 
which is rather consistent with an NEI (under-ionized) 
picture than an over-ionized one. 

For Reg 2, we examined the spectrum with a CIE model 
affected by the absorption of ISM. This model reproduced 
the spectrum with a x^/d.o.f of 200/160. The obtained 
spectrum and parameters are shown in figure 3c and ta- 
ble 2. We also tried an NEI model. The results were 
significantly improved to a x^/d.o.f of 168/159, and the 
normalization of the Gaussian at 1.2 keV became 0. The 
best-fit spectra and the parameters arc shown in figure 3d 
and table 2. We also checked the probability of over- 
ionization with the same procedure as used for Reg 1. 
As a result, T^/T^ was obtained to be 0.79lg'|!)g and then 
the plasma of Reg 2 is not over-ionized. Thus, we again 
adopted the NEI model as the best-fit for Reg 2. The 
plasma properties are essentially the same as Reg 1. 

The obtained results of Reg 1 and Reg 2 are incon- 
sistent with that obtained by Chandra (Koo et al. 2005). 
They reported that the spectra were represented by a CIE 
plasma with temperature of 0.3-0.5 keV and enhanced Zo,. 
This might be caused by the inappropriate background es- 
timation without considering the latitudinal dependence 
of GRXE in the Chandra data analysis. 

3.3. Hard X-Ray Emission 

3.3.1. Background Modeling 

The hard X-ray emission widely extends around molec- 
ular clouds associated with W51B but the distance to the 
emitting region is unclear. We can estimate the distance 
from iVn. Thus, the accurate estimation of the back- 
ground emission toward Reg 3 is required. We first mod- 
eled the spectrum of the BCD region and then modified 
it to Reg 3, considering the difference of the Galactic ab- 
sorption column densities and the intensities of GRXE. 
The modeling method is described below. 

As mentioned in subsection 3.2.1, the emission toward 
the BCD region consists of NXB, CXB, and GRXE. 
Unlike subsection 3.2.1, we here modeled the GRXE spec- 
trum by a three temperature CIE model (Uchiyama et al. 
2009). A Gaussian was added at 1.2 keV to compensate 
the lack of the Fe L-shell lines (see subsection 3.2.2). We 
also considered the Fe I Ka line of GRXE (Ebisawa et al. 
2008) to represent the emission in the hard X-ray band 
more accurately. Because the lowest temperature compo- 
nent of GRXE is considered as a foreground emission, the 
absorption for it was independent of that for the other 
two CIE components. After subtracting NXB. we fitted 
the spectrum of BCD, as shown in figure 4, using the fol- 
lowing spectral model. 
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Fig. 4. XIS spectrum of the BGD region with the best-fit 
models in equation (1). For visibility, data and models for 
only the FI are displayed. Black data shows the XIS data 
while a brown curve shows the CXB component. Magenta, 
blue, cyan, green, and purple curves show the best-fit mod- 
els of CIEl, CIE2, CIE3, Gausl and Gaus2 components in 
equation (1), respectively. 

Absml X Powerlawl -t- Absm2xCIEl 

+ Absm3 X (CIE2 + CIE3 + Gausl + Gaus2). (1) 

Powerlawl, CIEl, CIE2, and CIE3 represent CXB, FE, 
the medium and high-temperature component of GRXE. 
Gausl and Gaus2 are a Gaussian at 1.2 keV and 6.4 kcV 
lines, respectively. Absml, Absm2, and Absm3 represent 
the ISM absorptions for CXB, FE, and medium-f high- 
temperature components of GRXE, respectively. The 
spectral parameters for CXB including Absml were fixed 
at those used in subsection 3.2.1. GRXE has an axial 
symmetric structure for the Galactic center and its lon- 
gitudinal scale length of 45° (Uchiyama et al. 2013) is 
smaller than that of absorbing gas of ~100° (Amores & 
Lepine 2005). Therefore, we assumed that the GRXE 
emitting plasma is located at half distance to the edge of 
the Galaxy along the line of sight and the gas is divided 
to near and far side components for the plasma. Hence, 
we set Absm3 = 1/2 x Absml. For FE, we treated Absm2 
as a free parameter. The elemental abundances of CIEl 
were fixed at the solar value. Z-pe of CIE2 and CIE3 were 
set to be common and the other abundances were fixed 
at the solar value. The fitting model with equation (1) 
well represented the spectrum with a x^/d.o.f of 186/167 
(sec figure 4), where the best-fit parameters are listed in 
table 3. 

To derive the background emission in Reg 3, we scaled 
the normalization of CIEs and Gaussians of the BGD 
region by a factor of 1.93, which was obtained by the 
same estimation as used in subsection 3.2.1. For Absml 
and Absm3, we adopted the column densities of through- 
out the Galaxy and half of it toward Reg 3, respectively 
(3.3x10^2 cm-2 and 1. 7x10^2 cm'^). The TVh of Absm2 
and powerlaw model for CXB were assumed to be the 
same as those of BGD. Then, we obtained the background 
model for Reg 3. 



No. 
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Table 3. The best-fit result of the BGD region. 



Parameter 



Foreground 
(CIEl) 



Nh (10^^ cm-^) 0.49 
Te (keV) 
Zpe (solar) 

0.5- 



Mcdium 
(CIE2) 



-10 



/( 



6.4 



-0.12 

11 +°°^ 
"■^^ -0.03 

1.00 (fixed) 
1.60 (52.9) 



0.58 (fixed) 
49 +0 09 

"■^^ -0.13 



High 
(CIE3) 



o Qo +1.49 
■"••"^ -0.86 

^•"J^ -0.19 



2.25 (11.7) 7.25 (10.6) 



6.16 



6.16 



XVd.O.f 



186/167 
~ke^ 



Observed flux in tEe 0.5-10 keV band Hi units of 
10~^^ erg s~^ cm~^. Values in parentheses are the absorption- 
corrected values. 

^ Observed photon flux of 6.4 keV line in units of 10~^ pho- 
tons s"'^ cm~^. 



1 1 1 1 1 




' ' 1 


_ 1 1 1 1 1 







0.5 1 2 5 10 

Energy (keV) 



Fig. 5. XIS spectrum of Reg 3 with the best-fit model of 
equation (2). For visibility, data and models for only the 
FI are displayed. Black data shows the XIS data. Red, or- 
ange, green, and blue curves show the models for Background 
model, NEI, Gaus3 and CIE4 components in equation (2), 
respectively. 

3.3.2. Reg 3 

To investigate the nature of the extended hard X-ray 
emission, we analyzed Reg 3 as shown in figure 2b. The 
bright hard X-ray sources, the PWN candidate, and two 
Hll regions G49.0-0.3 and G49.2-0.3 were excluded from 
the region. 

We show the NXB-subtracted spectrum of Reg 3 in fig- 
ure 5. In this region, the soft thermal X-rays from W51C 
also contribute as background emission. We represented 
the W51C emission by the NEI model of Reg 1 (see sub- 
section 3.2.2). We did not use the Reg 2 results, since the 
errors are larger but the results are consistent with those 
of Reg 1. We tried to fit the hard X-ray emission with a 
CIE or a powerlaw model as follows. 

Background model + Absm4x (NEI + Gaus3) 

+ Absm5xCIE4, (2) 

or 

Background model + Absm4 x (NEI + Gaus3) 



+ Absm5 X Powerlaw2. (3) 

Here, the background model consists of CXB and the 
scaled GRXE. NEI and Gaus3 are the W51C emission 
and the 1.2 keV line making up for the lack of Fe-L shell 
lines, respectively. We set the normalization of the NEI 
and Gaus3 to be free. The iVn of Absm4 and Absm5, the 
normalization. To, and Zpo of CIE4, and the normalization 
and photon index of Power law2 were also free parameters. 
The fitted spectrum and parameters are shown in figure 5 
and table 4. 

The obtained value of A^h of Absm5 is smaller than 
~3xl0^^ cm~^, the total column density toward this di- 
rection, suggesting that the hard X-rays is Galactic ori- 
gin. This value is similar to ~1.5xl0^^ cm~^ derived by 
the X-ray spectral analysis of Hll regions, G48.9— 0.3 and 
G49. 0—0.3, located at the front side of the clouds (Koo 
et al. 2002). Therefore, we consider that the hard X-ray 
emitting region is located around the clouds and W51C, 
although this is not conclusive due to the systematic un- 
certainties stemming from the estimation of GRXE and 
spectral modeling of SNR emission. 

The analyzed spectrum of Reg 3 contains emission 
from point sources associated with W51B and background 
sources. We estimated the contribution of these point 
sources to the extended emission as described below. 

First, we investigated a time variability of the X-ray 
emission from point sources. The image of FIs in the 
2.5-8.0 keV band was compared with the Chandra im- 
age convolved with the PSF of the Suzaku XRT. We 
found a compact source, whose flux is similar to that of 
G49. 0—0.3, in the convolved Chandra image, located at 
{I, b) = (49.°0206, -0.°2629), while there is apparently no 
object in the FIs. It might be a compact source such 
as an active galactic nucleus like a blazar or a flare star 
which have a short term variability. We did not find other 
differences between the two images. 

Next, we estimated the total X-ray luminosity of the 
point sources in Reg 3 using the Chandra data. The to- 
tal photon counts of the point sources in Reg 3 without 
the above bright source were converted to fluxes in the 
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Table 4. The fit results of Reg 3* 



Parameter 



NEIT 
Nh (10 



Model of equation (2) Model of equation (3) 



22 



VEM (10^^ cni~^)^ 



2.42 (2.41) tVio 
5.02 (5.05) t'oil 



2.40 (2.38) 1^12 
4.66 (4.86) 



Thermal (CIE4) 
Nn (10^2 

VEM (10^8 cm-3)t 
To (keV) 
Zpc (solar) 
Luminosity^ 



0.86 (0.98) toll 
6.79 (6.50) 
5.08 (4.93) 
0.24 (0.23) 
6.53 (6.14) 



-0.86 
'0.75 
+ 1.12 
-0.83 
+0.16 
-0.16 
+0.83 
-0.72 



Nonthermal (Powerlaw2) 
iVn (1022 cm-2) 
Photon index 
Luminosity^ 



1.15 (1.30) 



+0.25 
-0.19 

2.23 (2.29) toll 
7.02 (6.63) tlfe 



xVd.o.f 443/375 445/376 

* 'Ihc values in the parentheses arc obtained by increasing the GKXE spec- 
trum by a factor of 1.25, for systcmatics. 

''' The parameters except for VEM are fixed at the values of the NEI model 
for Reg 1. 

Volume emission measure, VEM = J ncn^idV, where ric and tih arc the 
electron and Hydrogen densities, respectively, and V is the emitting volume, 
assuming the distance of 6 kpc. 



§ Unabsorbed luminosity in the 2.0-10.0 keV band in units of 10'^'^ 
for the distance of 6 kpc. 



erg s" 



2.0-10 keV band by using WebPIMMS2 assuming the emis- 
sion from a CIE plasma attenuated by absorption of ISM. 
We assumed 3 keV and the solar values for the tempera- 
ture and abundances of the CIE model, respectively, con- 
sidering an emission model for typical OB-type stars and 
young stellar objects (Schulz et al. 2001). Because the dis- 
tances to the sources are unknown, we adopted the A^h of 
2.1xl022 cm~2^ which is an intermediate value obtained 
by the X-ray analysis toward W51B (Koo et al. 2002). 
Using the obtained unabsorbed flux, the net luminosities 
of the point sources are 1.02xl0'^"^ erg s~^. 

We cannot distinguish which model (thermal or non- 
thermal) is statistically plausible for the present data. The 
total flux of the point sources toward Reg 3 has only minor 
contribution to that of the hard X-rays using the back- 
ground model, even considering the systematic error for 
the intensity profile of GRXE (see subsection 3.2.1). 

4. Discussion 

Contribution of Suprathermal Electrons to 
Ionization of Thermal Plasma 

In SNRs, electrons injected into the first order Fermi 
acceleration are expected to have suprathermal energies. 
Because such electrons can effectively ionize ions in ther- 
mal plasma, we can consider the possibility that the elec- 
trons make the plasma over-ionized (e.g., Kato & Masai 
1992). Cosmic rays responsible for gamma-ray emission 
in SNRs are accelerated around the shock whereas X-ray- 



emitting plasmas are formed in the downstream of the 
shock. Thus the cosmic rays and the plasmas would in- 
teract with each other. Note that we cannot resolve two 
regions due to the lack of angular resolution of a cur- 
rent gamma-ray telescope. The GeV emitting region in 
W51C overlaps with Reg 1 and Reg 2 (Abdo et al. 2009). 
However, we could not find any significant features of the 
over-ionized plasma in this SNR (subsection 3.2.2). One 
possibility is that the suprathermal electrons exist but 
the elapsed time of the ionization is too short to make 
the plasma over-ionized. To examine this idea, we here 
estimate the elapsed and characteristic timescales of ion- 
ization by suprathermal electrons. 

We assume the energy distribution of the suprathermal 
electrons to be quasi-thermal, which means it is roughly 
a Maxwellian. We also assume the "temperature" of the 
electrons to be 10 keV. In this case, the characteristic 
timescale of ionization by suprathermal electrons can be 
represented by (crv)"^, where a and v are the cross sec- 
tion of collisional ionization and the velocity of suprather- 
mal electrons, respectively. We consider the ionization of 
H-like Mg ions. By retrieving the value of a from the 
ALADDIN database'^, {ctv)~^ for the above situation is de- 
rived to be ~5xl0^° cm~^ s. 

The ionization parameter of suprathermal electrons is 
n'J^' , where rig and t' are the number density of suprather- 
mal electrons and the elapsed time of the ionization. 
Because n'^ is unknown, we assume that n'^ scales ng, the 
number density of thermal electrons. We take the den- 



littp: / /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 



(http: / /dpc. nifs.ac.jp/aladdin) 



No. ] 

sity ratio, rig/rie, to be 10"^ obtained around the shock 
front of SN 1006 where efficient acceleration occurs. The 
ratio should be taken as the upper limit for the middle- 
aged SNR W51C since Bamba et al. (2003) derived the 
amount of nonthermal electrons including suprathermal 
ones which have the energy of 10-100 keV. 

The value of the ionization parameter rict for the ther- 
mal electrons is about 2x10^^ cm~^ s from the spectral 
analysis (subsection 3.2.2). Then, by scaling this value 
for suprathermal electrons with the above n'^/ric, we de- 
rive n'^t' to be 2x10* cm~^ s under the assumption of 
t' ^ t, where t is the elapsed time of ionization by the 
thermal ones. The value of n[,t' is much shorter than the 
characteristic timescale of the ionization by suprathermal 
electrons. If the suprathermals are more abundant due to 
the second order Fermi acceleration processes (e.g., Bykov 
et al. 2000), n'^/n^ is about a few% (Masai et al. 2002). 
Even in this case, the value of n'^t' is still by one order of 
magnitude shorter than the characteristic one. Hence, the 
ionization by suprathermal electrons is negligible and the 
non-detection of an over- ionized plasma in W51C do not 
rule out the possibility that the cosmic-ray acceleration is 
ongoing there. Note that for a more realistic discussion, 
the time evolution of tt-o should be taken into account. 

4-2. Origin of Hard X-Ray Emission 

The hard X-ray emission in Reg 3 is represented by the 
thermal plasma model with ^^5 keV or the nonthermal 
emission with the photon index of ^2.2. The obtained A^h 
is quite smaller than that throughout the Galaxy, suggest- 
ing a Galactic origin. The emission probably have diffuse 
nature because the luminosity of 1x10'^^ erg s^^ in the 
0.5-10 keV band cannot be explained by the ensemble of 
the point sources located in the region. A PWN is un- 
likely to be the origin because no pulsars have been found 
in Reg 3. Based on these results, we discuss the possible 
origins of the diffuse hard X-ray emission. 
4.2.1. Stellar Winds from OB stars in W51B 

The hard X-ray emission overlaps the molecular clouds 
associated with the star-forming region W51B. There are 
at least four Hll regions around Reg 3. Such an envi- 
ronment leads us to the idea that the hard X-rays orig- 
inate from shock-induced bubbles due to stellar winds 
from OB stars (e.g.. Weaver et al. 1977). Stellar wind 
bubbles would have thermal and/or nonthermal emission 
(e.g., Wolk et al. 2002; Ezoe et al. 2006). 

Wc first consider the possibility that the hard X-ray 
photons originate from shock-heated thermal plasma via 
the stellar winds. The volume of the emission region is 
approximated to be a sphere of 10-pc radius estimated 
from the apparent size of Reg 3 and the line of sight 
depth of the clouds (Garpenter & Sanders 1998). From 
the equation (3) of Ezoe et al. (2006), the total thermal 
energy with the temperature of 5 keV is obtained to be 
■~2x 10^" erg. As the age of the H II regions around Reg 3 
is about 1 Myr (Kim et al. 2007), the required power is 
--6x10^'^ erg s~^ According to Kim et al. (2007), the 
spectral types of the most massive stars in H II regions 
G49.0-0.3 and G49.2-0.3 are 09 and 04, respectively. 
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Muijres et al. (2012) estimated that the mass loss rate and 
the velocity of the stellar winds of an 04 star is about 
10~5-^ M0 yr~^ and 4000 km s~^, respectively. Then 
the kinetic power of the stellar winds is calculated to be 
^10^'' erg s~^ from the equation (6) of Ezoe et al. (2006). 
Therefore, thermal emission from stellar winds might be 
the origin of the hard X-rays. 

Next, we discuss nonthermal origin with stellar winds; 
bremsstrahlung from nonthermal (> 10 kcV) electrons, in- 
verse Gompton (IC) scattering by hundreds-of-MeV elec- 
trons, and synchrotron emission from TeV electrons. We 
can rule out the bremsstrahlung emission on the basis of 
the spectral index. In a dense environment (the gas den- 
sity > 10 cm~^), bremsstrahlung emission dominates the 
hard X-ray band. The spectral index of energetic elec- 
trons should be flatter than the initial spectrum because 
lower energy electrons lose their energy more quickly than 
more energetic ones via the Goulomb loss. The resulting 
X-ray emission has a spectral index of ~ 1.2 (Uchiyama 
et al. 2002), but the observed index is '^2.2. 

The possibility of IC scattering can also be ruled out. 
Electrons with energies of several hundred MeV (with a 
Lorentz factor of lO^-lO'^) can produce hard X-ray emis- 
sion by scattering the IR emission from star lights in 
W51B. The steep hard X-ray spectrum could be the result 
of synchrotron loss. However, the cooling time of the elec- 
trons even with the magnetic field of 100 /iG is ^10 Myr, 
which is much longer than the age of W51B. 

Synchrotron emission from tens of TeV electrons pro- 
duced by stellar winds could explain the diffuse X-rays. 
We roughly estimated the amount of the energy of the 
TeV electrons responsible for the X-ray emission in the 
0.5-10 kcV band assuming a magnetic field of 10 /iG and 
then it is obtained to be about 1 x 10**^ erg. Since the cool- 
ing time of the TeV electrons is a few 10'^ yr, the energy 
should be injected within the timescale. Then, the aver- 
age energy input is constrained to be larger than about 
10'^''^ erg s~^, which can be supplied by stellar winds from 
the OB stars around W51B. 
4.2.2. SNR W51C 

The hard diffuse X-ray emission also overlaps the west- 
ern part of the shell of W51G. Thermal emission from the 
SNR is less likely to be the origin of the diffuse X-rays 
since the typical plasma temperature of SNRs is below a 
few keV. Therefore, nonthermal origin is more preferable 
interpretation. We can immediately rule out the possibil- 
ity of bremsstrahlung and IG scattering because the spec- 
tral index of the X-rays is soft (~2.2: subsection 4.2.1). 

If particle acceleration terminates at the early Sedov 
stage, synchrotron emission from tens of TeV electrons 
cannot explain the hard X-rays since the cooling time 
of the high energy electrons is a few 10'^ yr under a 
magnetic field of '^10 fjiG. Alternatively, the observed X- 
ray flux might be explained by the synchrotron emission 
from secondary particles produced by decay of charged pi- 
ous generated by interaction of cosmic rays in dense gas. 
However, the synchrotron X-ray emission from secondary 
particles is fairly dim in the hard X-ray band even with 
higher gas density than 100 cm^"^ in an SNR having a 
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similar 7-ray spectrum and an age to W51C (e.g. Ajello 
et al. 2012). 

On the other hand, if particle acceleration takes place 
over the SNR age, tens-of-TeV electrons can exist in 
middle-aged SNRs (Sturner et al. 1997; Nakamura et al. 
2012) and then synchrotron emission from the TeV elec- 
trons can explain the X-ray spectrum. The detection of 
possible synchrotron X-rays in the middle-aged SNR W44 
might be related to TeV electrons (Uchida et al. 2013) 
and supports the above scenario. Assuming the acceler- 
ation over the SNR age, very efhcient acceleration is un- 
likely to operate because n'^t' exceeds the characteristic 
timescale of the ionization by suprathermal electrons un- 
der a efficient acceleration (see subsection 4.1). For the 
acceleration only in the late phase of the SNR, the particle 
acceleration mechanism might be related to the interac- 
tion of the SNR shells with molecular clouds. Inoue et al. 
(2009) demonstrated that shock-compressed shells become 
turbulent owing to the inhomogeneity of preshock density 
and then amplification of magnetic field takes place in 
the shell. Particles could be accelerated at the secondary 
shocks arising from the turbulent flows. However, the sec- 
ondary shocks with a velocity of hundreds km s~^ cannot 
accelerate electrons up to TeV energies due to fast syn- 
chrotron cooling. 

An alternative idea is that high energy electrons are 
accelerated by SNR(s) in tenuous ISM, where the SNRs 
would effectively accelerate particles without deceleration 
of the shocks (e.g.. Tang & Wang 2005). Such an environ- 
ment might be provided by core-collapse supernova explo- 
sions before the birth of W51C. This idea is supported by 
the enhanced Mg abundance in Reg 1 and Reg 2. Because 
the main component of the plasma is possibly ISM (Koo 
et al. 2002 and this work; see Zsi and Zg in table 2), 
the enhanced Mg would come from the past explosion of 
20-25 M© stars, which supply a large amount of Ne and 
Mg (Tsujimoto et al. 1995). The existence of two PWN 
candidates around W51C provides another clue for the 
past core-collapse supernovae in this region. 

5. Conclusion 

We investigated the X-ray emission toward the middle- 
aged SNR W51C and its environment with XIS onboard 
Suzaku. The soft X-ray emission from the plasma of 
W51C is well represented by an optically thin thermal 
plasma in NEI with the temperature of ^0.7 keV with 
the enhanced Mg abundance. Although the thermal emis- 
sion region overlaps with the GeV emitting region where 
cosmic rays are accelerated in SNRs, there is no hint of 
an over-ionized plasma. These results are consistent with 
the situation that the elapsed time of the ionization by 
suprathermal electrons is too short to make the plasma 
over-ionized. 

The hard X-ray emission is spatially coincident with 
the molecular clouds associated with W51B. The spec- 
trum of Reg 3 is represented by an optically thin ther- 
mal plasma with a temperature of '^5 keV or a powcrlaw 
model with a photon index of ~2.2. The emission have 



diffuse nature because the luminosity of 1 x 10 erg s~^ in 
the 0.5-10 keV band cannot be explained by the ensemble 
of point sources. The smaller iVn than that throughout 
the Galaxy suggests a Galactic origin but a PWN is un- 
likely the candidate because no pulsars have been found 
in Reg 3. The origin of the hard X-rays can be interpreted 
with the stellar wind scenarios. The X-ray emission can 
be also interpreted by the synchrotron X-rays due to the 
TeV electrons continuously accelerated by the SNR W51C 
after the free expansion phase or in tenuous ISM while it 
cannot be explained by the synchrotron emission either 
of leptonic or hadronic process of particles accelerated by 
the early Sedov phase. 
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